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In  this  paper,  thermal  effect  of  a  diode-end-pumping  Tm:LuAG  crystal  is  analyzed  and  the  lasing 
characteristics  at  room  temperature  are  also  investigated.  Based  on  a  thermal  model,  the  thermal  focal 
length  of  Tm:LuAG  crystal  is  calculated  at  different  pump  energies  between  50  mj  and  300  mj.  The 
sensitivity  factor  of  the  thermal  lens  M=d(L/f)/dPabs  is  estimated  to  be  0.36m-1/VV  or  36m_1/J-  The 
thermal  focal  length  of  Tm:LuAG  crystal  is  also  measured  at  different  pump  energies.  Results  show  that 
the  experimental  data  are  very  well  in  accordance  with  the  theoretical  calculations.  Laser  characteristics 
of  Tm:LuAG  crystal  at  room  temperature  are  also  investigated.  Under  a  pump  energy  of  300  mj,  the 
maximum  output  energy  of  41.0  mj  is  obtained,  with  a  slope  efficiency  of  20.8%.  The  build-up  time  of  the 
pulse  is  2.76  ms,  and  the  beam  quality  factor  is  measured  to  be  M2=1.17. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

All-solid-state  lasers  operating  at  2  pm  are  useful  for  the 
applications  in  coherent  Doppler  LIDAR,  differential  absorption 
LIDAR  and  so  on  [1-4].  Thulium-doped  Lu3A150i2  (TrmLuAG) 
crystal  is  one  of  the  attractive  materials,  which  is  isomorphic  to 
YAG  [5].  It  has  not  only  the  advantages  of  higher  heat  conductivity, 
but  also  the  lower  population  density  at  the  lower  laser  level.  The 
emission  wavelength  of  TrmLuAG  (2.023  pm)  is  more  closer  to  the 
optical  transmission  window  than  that  of  Tm:YAG  (2.015  pm). 

In  1995,  diode-pumped  TrmLuAG  laser  with  an  optical-to-optical 
conversion  efficiency  of  23.6%  was  reported  [6].  In  2000,  an  output  of 
0.92  W  diode-pumped  TrmLuAG  laser  was  obtained  [7].  In  2004,  a 
single  mode  TrmLuAG  laser  was  developed  [8  .  In  2008,  the  max¬ 
imum  output  power  of  4.91  W  TrmLuAG  laser  with  a  slope  efficiency 
of  25.39%  was  obtained  [9].  In  2009,  a  single  frequency  TrmLuAG 
laser  was  reported  with  output  power  of  148  mW  [10].  Although  the 
laser  characteristics  of  TrmLuAG  laser  had  been  demonstrated,  few 
research  works  were  carried  out  on  the  thermal  effect  analysis  of  Tm: 
LuAG  crystal.  As  a  quasi-three-level  laser  system,  the  thermal  effect 
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of  Tm3+ -doped  crystal  could  not  be  ignored.  In  2011,  Yao  et  al.  [11] 
reported  the  thermal  analysis  of  TrmYLF  crystal  and  TrmYAP  crystal. 
In  2012,  Gaponenko  et  al.  [12]  reported  the  thermal  lensing  of  Tm3+- 
doped  KY(W04)2  crystal.  The  results  showed  that  the  thermal  effect 
of  Tm3+ -doped  crystal  was  quite  severe,  which  would  affect  the 
output  laser  characteristics. 

This  letter  focuses  on  the  analysis  of  thermal  effect  about  Tm: 
LuAG  crystal  and  the  lasing  characteristics  of  TrmLuAG  at  room 
temperature.  The  thermal  model  of  TrmLuAG  crystal  was  built  up. 
And  the  measured  thermal  focal  length  was  nearly  the  same  as  the 
calculated  value.  For  example,  the  thermal  focal  length  of  Tm: 
LuAG  crystal  was  92.8  mm  theoretically  and  102.7  mm  experi¬ 
mentally  under  the  same  conditions  of  pumping  source  and  the 
same  parameters  of  laser  crystal.  Under  pump  energy  of  300  mj, 
maximum  output  energy  of  41.0  mj  was  obtained,  with  slope 
efficiency  of  20.8%.  The  build-up  time  of  the  pulse  was  2.76  ms. 
And  the  beam  quality  factor  was  M2  =  1.17. 


2.  Thermal  effect  analysis  of  TmiLuAG  crystal 

2.1.  Thermal  model 

A  diagrammatic  sketch  of  diode-end-pumped  TrmLuAG  rod  is 
shown  in  Fig.  1.  This  model  was  similar  to  our  previous  report 
[13,14].  In  Fig.  1,  r  and  z  are  the  radial  and  the  axial  coordinates, 
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Fig.  1.  Model  of  pulsed  laser-diode  end-pumped  a  Tm:LuAG  rod. 


respectively,  b  is  the  radius  of  rod,  and  L  is  the  length  of  rod.  Pump 
light  is  propagated  in  the  z  direction. 

Along  the  z  direction,  the  pump-light-intensity  is  described  by, 
2  p 

I(r,z,t)  =  — Texp(-2r2/co2)exp[-fiz]G(t),  0 <z<L  (1) 

7TGJZ 

here,  P  is  the  peak  power  of  a  pump  pulse,  go  is  the  radius  of  the 
pump  beam  waist,  and  /?  is  the  absorption  coefficient  of  TrmLuAG 
rod.  G(t)  is  given  by, 

G(t)  =  [square(2^/t,  r  •  /)  + 1  ]/2  (2) 

square  is  the  rectangular  square  wave  function,  t  is  the  time  and 
/  is  the  pump  repetition  rate. 

The  heat  source  function  is  expressed  as, 

qv(r, z,  t)  =  fhjHr, z,  t )  =  2^-exp(  -2 r2 /co2)exp[-fiz]C(t)  (3) 


here,  rj  is  the  heat  fraction. 

In  single  pulse  duration  time,  it  satisfies  the  transient  heat 
conduction  mode  as  shown  in  Eq.  (4). 
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In  the  following,  an  integral  transform  method  is  used  to  get 
the  expressions  of  Ti  and  T2. 
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The  unknown  variables  are  described  as  follows, 
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The  boundary  condition  is  shown  in  Eqs.  (5)-(7). 
T\  =  Tw,r=b 

-k— +hT,  =hTe,z  =  0 

dZ 


k—+hTi=hTe,z  =  L 
dz 

Initial  condition  is  shown  in  Eq.  (8). 

Ti  =T0,  t  =  0,  0<z<L0<z<L,  0<r<b 

In  single  pulse  interval  time,  it  satisfies  the  transient  heat 
conduction  mode  as  shown  in  Eq.  (9). 
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,  r<t<T,  0  <  z  <  L,  0<r<b  (9) 


Boundary  condition  is  expressed  by  Eqs.  (10)-(12). 
T2  =  Tw,  r  =  b 


(10) 


R0(fJm,r)=J0(/3mr)  (20) 

_ 1 _ = _ - _  (21) 

Wm)  b2J'2(fimb)  (  j 

here,  J0(fimr)  is  the  zero  order  Bessel  function.  Eigenvalue  pm  is  the 
positive  root  of  the  equation  J0(JJmb)  =  0. 

Z(rjp,z)  =  rjp  cos  rjpZ+(h/k)  sin  rjpz  (22) 


[t]j  +  (h/k)2] 


L+ 


(h/k) 
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+  (h/k) 


(23) 


here,  rjp  is  the  positive  root  of  the  equation  tan  rjpL  =  (rjp(h/k+ 
h/k)/r]2-(hh/kk)). 

After  n  (n  >  2)  pump  pulses,  the  expression  of  temperature 
distribution  in  repeated  pulse  duration  time  could  be  given  as 
follows: 


-l<d^+hT2=hTe,  z  =  0  (11) 

dZ 

\i—+hT2  =  hTe,z  =  L  (12) 

dZ 

The  initial  condition  is  given  in  the  Eq.  (13). 

T2(r,z,t)\t  =  r  =  Ti(r,z,r)  (13) 

here  p  is  the  crystal  density,  C  is  the  specific  heat  capacity,  h  is  the 
heat  transfer  coefficient  of  the  air,  k  is  the  heat  conductivity.  Tw  is 
the  temperature  of  cooling  water.  Te  is  the  initial  temperature  in 
the  crystal.  T0  is  the  room  temperature. 
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After  it  is  pumped  by  n(n  >  2)  pulses,  the  temperature  distribu¬ 
tion  in  repeated  pulse  interval  time  is  given  by 


Jr  =0 


f,  r'Ro(Pm,r')Z(rip,z')- 
Jz  =  0 


OO  O op- 
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Ro{Pm,r)Z(r]v,z)- 
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eay£  ^  f)dt 


(25) 


In  the  rod,  the  thermally  induced  optical  path  difference  (OPD) 
is  induced  by  the  change  of  refractive  index  with  temperature  [15]. 
The  OPD  in  the  slice  z  to  dz  is  given  by 


dOPD(r,z,  t)  = 


^+(n-l)(w+l)aT 


:  [T (0, z,t)-T (r, z, t)]dz  (26) 


f  = 


The  thermal  focal  length  varied  with  time  is  expressed  as 

r2  r2 

1  e  _ ‘_e _ 

20PD  2  { [af + (”  - 1  + 1  )«t]  x  faL[T(0,z,t)-T(rb,z,t)]dz} 

here  re  is  the  effective  radius  of  pump  light. 


(27) 


2.2.  Simulation  results  and  analysis 

The  parameters  used  in  simulation  is  shown  in  Table  1. 

The  transient  temperature  distribution  and  the  time-varying 
thermal  focal  length  in  TrmLuAG  crystal  are  shown  in  Figs.  2  and  3  . 

As  shown  in  Fig.  2,  the  temperature  distribution  appeared 
jagged  with  time.  The  maximum  temperature  of  TrmLuAG  crystal 
arrived  at  72.7°C.  And  the  interval  that  axially  transient  tempera¬ 
ture  distribution  reached  steady  state  is  1.2  s.  In  Fig.  3,  the  time- 
varying  thermal  focal  length  appeared  jagged  and  decreased  with 
time.  The  thermal  focal  length  is  in  the  range  of  92.8-381.0  mm.  At 
pump  peak  power  of  30  W,  pump  frequency  of  50  Hz  and  pump 
pulse  duration  of  10  ms,  the  thermal  focal  length  is  381.0  mm  at 
the  beginning  and  stable  at  92.8  mm  after  about  120  pumped- 
cycles.  From  Figs.  2  and  3,  it  was  obvious  that  although  the 
pumping  source  was  a  pulsed  laser  diode,  the  temperature 
distribution  and  focal  lens  of  thermal  lensing  stabilized  to  a  steady 
state  while  pump  time  continued. 


Table  1 

Parameters  used  in  simulation. 


Parameters 

Value 

Temperature  of  cooling  liquid/K 

291 

Initial  temperature/K 

295 

Heat  transfer  coefficient  of  air/(W  cm-2  K-1  ] 

) 

8  x  10-4 

Heat  conductivity  of  Tm:LuAG/(W  cm-1  K-1) 

[16] 

0.079 

Pump  beam  waist/cm 

0.12 

Pump  peak  power/W 

30 

Pump  frequency/Hz 

50 

Pump  pulse  duration/ms 

10 

Heat  fraction 

0.6 

Crystal  radius/mm 

1.5 

The  length  of  TrmLuAG  rod/mm 

8 

Density  of  Tm:LuAG/(g/cm3)  [16] 

6.32 

Specific  heat  capacity  of  Tm:LuAG/(J  Kg-1  K-1 

)  (similar  toTrmYAG) 

590 

Absorption  coefficient  of  4  at%  Tm:LuAG/cm“ 

i 

3 

Thermo-optical  coefficient  (dn/dT)/(K_1)  [16] 

8.9  x  10-6 

Refractive  index  of  TrmLuAG  rod 

1.81 

Poisson's  ratio 

0.3 

Thermal  expansion  coefficient/(I<_1)  [16] 

7.9  x  10“6 

100 

80 


40 

20 


Fig.  2.  Axially  transient  temperature  distribution  in  TrmLuAG  crystal. 


Fig.  4  gave  the  thermal  focal  length  of  TrmLuAG  crystal  under 
different  pump  energies  (pump  peak  power  x  Pump  pulse  dura¬ 
tion).  And  the  sensitivity  factor  of  the  thermal  lens  M=d(\lf)ldPabs 
was  about  0.36  m -1/W  or  36  m_1/J-  The  smaller  the  pump  spot 
radius  was,  the  shorter  the  thermal  focal  length  became.  There¬ 
fore,  the  value  of  M  would  be  bigger,  which  meant  that  the  laser 
crystal  would  be  more  sensitive  to  the  thermal  effect. 
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3.  Experimental  setup 

Fig.  5  shows  the  experimental  setup  of  the  Tm:LuAG  laser.  The 
oscillator  was  a  bowtie  cavity  in  a  bidirectional  fashion.  The  pump 
source  was  pulsed  fiber-coupled  laser-diode  (LD)  with  a  peak 
power  of  30  W  and  central  wavelength  of  788  nm.  The  frequency 
and  the  pulse  duration  of  the  LD  was  50  Hz  and  10  ms,  respec¬ 
tively.  The  fiber  core  had  a  diameter  of  400  pm  and  a  numerical 
aperture  of  0.22.  The  output  beam  of  the  laser-diode  was  collected 
by  25  mm  focal  length  collimating  lens,  and  it  was  focused  into  the 
TrmLuAG  crystal  using  a  lens  with  75  mm  focal  length.  Then,  the 
radius  of  pumping  beam  waist  was  0.6  mm.  The  TrmLuAG  laser 
crystal  had  a  Tm-doped  concentration  of  4%,  with  the  dimensions 
of  03  x  8  mm3.  Both  sides  were  coated  antireflection  near  788  nm 
(. R  <  0.5%)  and  2.02  pm  ( R  <  0.5%).  The  temperature  of  the  crystal 
was  kept  at  291  K  with  a  TE-cooler.  The  resonator  was  designed  to 
be  insensitive  to  the  thermal  lens  of  TrmLuAG  crystal  by  consider¬ 
ing  the  thermal  focal  length,  as  calculated  above.  The  radius  of 
oscillating  beam  {oo)  would  be  changed  because  of  the  thermal 
effect  of  TrmLuAG  crystal.  And  the  stability  of  the  cavity  would 
turn  to  bad.  Therefore,  the  dco/df  was  introduced  as  an  important 
factor  to  characterize  whether  a  cavity  was  designed  successfully. 
The  doo/df  for  our  experiment  was  approximate  equal  to  0  unless 
the  thermal  lens  was  shorter  than  80  mm.  Ml  was  a  high- 
reflectivity  flat  mirror.  M2  was  a  flat  output  coupler  with  a 
transmission  of  3.5%  at  2.02  pm.  M3  and  M4  were  curved  mirrors 
with  5  m  radius  of  curvature  and  coated  with  high-reflectivity  at 
2.02  pm.  The  total  resonator  length  was  0.3  m. 


4.  Experimental  results  and  discussion 

The  output  energy  of  TrmLuAG  laser  was  measured  by  an 
energy  meter,  as  shown  in  Fig.  6.  The  threshold  of  the  laser  was 


95.3  mj.  The  maximum  output  energy  was  41.0  mj  with  pump 
energy  of  300  mj.  The  slope  efficiency  was  20.8%.  There  is  no 
saturation  phenomenon  of  TrmLuAG  output  energy.  It  meant  that 
the  cavity-design  was  successful. 

The  output  laser  was  detected  by  an  InGaAs  detector  connected 
with  a  digital  oscilloscope  (Tektronix  TDS3032B).  The  build-up 
time  of  the  pulse  was  2.76  ms  when  the  pump  energy  was  300  mj, 
as  shown  in  Fig.  7. 

At  the  maximum  output  energy,  a  beam  quality  factor  of 
M2  =  1.17  was  measured  by  the  traveling  knife-edge  method.  By 
the  way  of  the  lens-transformation,  as  shown  in  Fig.  8,  we  could 
know  the  position  of  the  beam  waist  and  the  beam  divergence 
angle  in  the  resonator,  according  to  the  following  equations:. 

«o  =  -  F  =0)0  (28) 

V0>-F)2+/2 


(7  =  F+ 


(f-F)F2 

(i-f)2+/2 


(29) 


where, /is  the  confocal  parameter,  /  =  jzcoq2 /L  l  or  l1  is  the  distance 
between  the  lens  and  the  waist  in  the  cavity  or  outside  the  cavity. 
F  is  the  focal  length  of  the  lens. 

The  position  and  the  radius  of  the  beam  waist  in  the  cavity 
would  be  different  by  considering  or  ignoring  the  thermal  lens  of 
the  TrmLuAG  crystal.  This  difference  should  disappear  by  changing 
the  value  of  the  thermal  focal  length.  However,  it  was  a  huge 
amount  of  calculation. 

Based  on  the  similar  principles,  by  means  of  software  about 
LasCAD,  the  thermal  focal  length  could  be  defined  easily.  Specifi¬ 
cally,  the  position  of  the  beam  waist  and  the  beam  divergence 
angle  in  the  cavity  could  be  determined  when  designing  the 
resonator.  Supposing  a  thermal  lens,  changing  the  value  of  the 
thermal  focal  length,  then  the  position  of  the  beam  waist  and  the 
beam  divergence  angle  in  the  cavity  would  be  changed 
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Fig.  6.  Output  energy  versus  pump  energy. 


Fig.  8.  Beam  waist  transformation  by  the  lens. 
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corresponding.  The  value  of  the  thermal  focal  length  changed  until 
the  position  of  the  beam  waist  and  the  beam  divergence  angle  in 
the  resonator  was  the  same  with  the  measuring  values.  The 
thermal  focal  length  determined  by  the  above  method  was  named 
as  experimental  results. 

The  thermal  focal  length  of  TrmLuAG  crystal  was  102.7  mm  at 
pump  energy  of  300  mj  (pump  peak  power  of  30  W  x  pulse 
duration  of  10  ms)  and  frequency  of  50  Hz.  By  this  method,  the 
value  of  thermal  focal  length  could  be  identified  at  pump  energy  of 
150  mj,  200  mj  and  250  mj.  Fig.  9  gave  the  calculated  results  and 
experimental  results  about  thermal  focal  length  versus  pump 
energy.  The  experimental  results  coincide  with  the  theoretical 
calculations.  The  deviation  might  due  to  the  missing  parameters  of 
TrmLuAG  crystal. 

5.  Summary 

In  summary,  we  reported  an  analysis  of  thermal  effect  about  Tm: 
LuAG  crystal  and  its  lasing  characteristics  at  room  temperature. 


Based  on  the  thermal  model  of  TrmLuAG  crystal,  the  thermal 
focal  length  was  estimated  theoretically.  By  means  of  software 
about  LasCAD,  the  thermal  focal  length  was  also  defined  from 
experimental  results.  The  experimental  results  coincide  with  the 
theoretical  calculations.  The  deviation  might  due  to  the  missing 
parameters  of  TrmLuAG  crystal.  Under  pump  energy  of  300  mj, 
maximum  output  energy  of  41.0  mj  was  obtained,  with  slope 
efficiency  of  20.8%.  The  build-up  time  of  the  pulse  was  2.76  ms. 
And  the  beam  quality  factor  was  M2  =  1.17. 


Acknowledgments 

This  work  is  supported  by  Science  and  Technology  Department 
of  Jilin  Province  in  China  (No.2012011113). 


References 


[1]  Y.F.  Li,  B.Q.  Yao,  Z.G.  Wang,  Y.Z.  Wang,  Y.L.  Ju,  Chin.  Opt.  Lett.  4  (2006)  470. 

[2]  B.Q.  Yao,  L.Q.  Dong,  Y.Z.  Wang,  Y.L.  Ju,  X.  Yu,  W.F.  Sun,  Y.F.  Li,  Acta  Opt.  Sin.  24 
(2004)  79. 

[3]  J.  Yu,  B.C.  Trieu,  E.A.  Modlin,  U.N.  Singh,  M.J.  Kavaya,  S.  Chen,  Y.  Bai,  P.J.  Petzar, 
Opt.  Lett.  31  (2006)  462. 

[4]  C.  Nagasawa,  T.  Suzuki,  H.  Nakajima,  H.  Hara,  K.  Mizutani,  Opt.  Commun.  200 
(2001)  315. 

[5]  J.D.  Kmetec,  T.S.  Kubo,  T.J.  Kane,  C.J.  Grund,  Opt.  Lett.  19  (1994)  186. 

[6]  N.P.  Barnes,  M.G.  Jani,  R.L.  Hutcheson,  Appl.  Opt.  34  (1995)  4290. 

[7]  V.  Wulfmeyer,  M  Randall,  A.  Brewer,  R.M.  Hardesty,  Opt.  Lett.  25  (2000)  1228. 

[8]  I<.  Scholle,  E.  Heumann,  G.  Huber,  Laser  Phys.  Lett.  1  (2004)  285. 

[9]  C.T.  Wu,  Y.L.  Ju,  Y.F.  Li,  Z.G.  Wang,  Y.Z.  Wang,  Chin.  Opt.  Lett.  6  (2008)  415. 

[10]  C.T.  Wu,  Y.L.  Ju,  Q.  Wang,  Z.G.  Wang,  F.  Chen,  R.L.  Zhou,  Y.Z.  Wang,  Laser  Phys. 
Lett.  6  (2009)  707. 

[11]  B.Q.  Yao,  P.B.  Meng,  G.  Li,  Y.L.  Ju,  Y.Z.  Wang,  J.  Opt.  Soc.  Am.  B  28  (2011)  1866. 

[12]  M.S.  Gaponenko,  P.A.  Loiko,  N.V.  Gusakova,  K.V.  Yumashev,  N.V.  Kuleshov, 
A.A.  Pavlyuk,  Appl.  Phys.  B  108  (2012)  603. 

[13]  G.Y.  Jin,  J.  Wu,  C.T.  Wu,  E.H.  Fan,  C.  Wang,  X.Y.  Chen,  Y.J.  Yu,  X.C.  Tan,  Chin. 
J.  Lasers  40  (2013)  (1002002-1). 

[14]  J.  Wu,  C.T.  Wu,  E.H.  Fan,  C.  Wang,  X.Y.  Chen,  Y.J.  Yu,  G.Y.  Jin,  Proc.  SPIE  8904 
(2013)  (890413-1). 

[15]  X.L.  Song,  B.B.  Li,  S.Y.  Wang,  D.F.  Cai,  J.G.  Wen,  Z.  Guo,  Chin.  J.  Lasers  34  (2007) 
1476. 

[16]  X.D.  Wang,  X.D.  Xu,  T.C.  Zang,  C.L.  Ma,  Z.W.  Zhao,  J.  Xu,  Chin.  J.  Rare  Earth  Soc. 
27  (2009)  745. 


